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O ; ABSTRACT 

I/"") , Large samples of Field Horizontal Branch (FHB) stars make excellent tracers of the 

■ Galactic halo; by studying their kinematics, one can infer important physical proper- 

ties of our Galaxy. Here we present the results of a medium-resolution spectroscopic 
survey of 530 FHB stars selected from the Hamburg/ESO survey. The stars have a 
mean distance of ~ 7 kpc and thus probe the inner parts of the Milky Way halo. We 
measure radial velocities from the spectra in order to test the model of Sommer-Larsen 
+-> ■ et al., who suggested that the velocity ellipsoid of the halo changes from radially- 

dominated orbits to tangentially-dominated orbits as one proceeds from the inner to 
the outer halo. We find that the present data are unable to discriminate between this 
model and a more simple isothermal ellipsoid; we suggest that additional observations 
towards the Galactic centre might help to differentiate them. 
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1 INTRODUCTION 

The stellar halo only comprises about 1 % of the luminous 
mass of the Galaxy and is composed of very old, metal- 
poor stars. In spite of this, the present-day dynamical and 
chemical state of the stellar halo has exerted considerable 
influence in shaping our understanding of the formation of 
large disk galaxies such as the Milky Way, particularly in its 
very early stages. 

Large samples of Field Horizontal Branch stars (FHBs) 
in the Galactic halo are e xcellent "test particles" for studies 
of halo kinematic s (e.g. ISommer -LarserietaLMlQQ^IlQQTl: 
ISirko et al.ll2004ah ; tracing the mass HlBe^r^^t^L||2u 04) and 
the merger history jKinman et aljlL993TBrowrTet^Lll2004T) 
of the Milky Way. Their intrinsic brightness and the relative 
ease with which distances can be obtained, together with 
their relatively clean spectra, also make them ideal probes of 
distances to High Velocity Clouds (e.g. ISchwarz et al.lll99<4 
Thom et al., in preparation). 

The HK survey of Beers and collaborators (e.g. 



iBeers et al.lll98,4 Il993l has identified ~ 12,000 candidate 
Horizontal Branch stars, of which about half are expected 
to be on the Horizontal Branch; the other half appear to 
be a mixt ure of A-type main-sequ ence stars and halo blue 
stragglers (Be ers et alJ fl988. 1996). More recently, several 
groups have reported large and clean samples of Horizontal 
Branch stars. The Hamburg/ESO survey (HES), originally 
designed to identify low-luminosity quasars, has been shown 
to contain a remarkable sample of interesting stars, includ- 
ing 8321 FHBs. This sample has a contamination level by 
non- FHBs (mostly high-gra vity A-type stars) of less than 
16% JChristlieb et al Jl2004h . A subset of these stars is used 
here to constrain models of halo kinematics. 

The Sloan Digital Sky Survey (SDSS), although pri- 
marily designed for extragalactic studies, has also revealed 
a large sample of 1 170 Horizontal Bran ch stars at distances 
up to ~ 100 kpc iSirko et al] 1200 4b'). From this sample 
ISirko et "all j2004al) measured an isotropic velocity ellipsoid 
for the outer halo. This result contrasts with halo stars in 
the solar neighbourhood, which show an ellipsoid elongated 
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Figure 1. Radial velocity differences, as measured using the 
H<5 and H7 lines, versus Signal-to-Noise ratio of the spectra, for 
our sample stars. The vertical and horizontal lines indicate the 
adopted S/N > 15 and l^diffl < lOOkms - 1 limits, respectively 
(see text). 
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Figure 2. RV1 versus RV2 for the 21 stars with repeated mea- 
surements. The solid line is the one-to-one correspondence. 



in the radial direction dChiba fc Beerd 120001 iGouldl I2003L 

Isooj. 

IWilhelm et al] il999T> have described a technique for 
separating Horizontal Branch stars from higher gravity A 
stars on the basis of broadba nd UBV colours and medium- 
resolution spectroscopy. IClewlev et al.lll2002l . l2004) explored 
similar techniques based on broadband colours and medium- 
resolution spectroscopy, and also presented a method rely- 
ing solely on medium resolution spectroscopy. They used 
this method to identify ~ 100 stars at distances of > 30 kpc 
with the aim of providing better constraints on the mass of 
the Milky Way. 

iBrown et ail i2004h used the full 2MASS point-source 
catalogue to select close to 100,000 objects, of which they 
expect 47 % to be on the Blue Horizontal Branch. They find 
an absence of structure in the inner halo (the sample has 
distance d < 9 kpc), concluding that there have been no 
major accretion events in the inner halo over the past few 
Gyrs. 

Here we present the results of a programme to ob- 
tain radial velocities (RVs), and study the kinematics, of 
more than 500 FHB stars selected from the HES cata- 
logue. We wish to test spe cifically the predictions of the 
ISommer -L^seii^^lJ |^94|) (h ereafte r, SLFC) model, as rc- 
fined bv ISommer-Larsen et al] lll997|) . In the SLFC model, 
the velocity ellipsoid of the halo changes from radially- 
dominated orbits to tangentially-dominated orbits as one 
proceeds from the inner to the outer halo. We attempt to 
distinguish between this model and the simpler isothermal 
halo model, which specifies an isotropic distribution of stel- 
lar orbits, independent of Galactocentric radius, in the outer 
halo. Section 2 presents the sample and observations. It also 
includes details on the radial velocity measurements and se- 
lection criteria. Section 3 discusses the models which we wish 
to test. Section 4 details the final sample selection, divi- 
sion into fields, and calculation of the line-of-sight velocity 
dispersions, ctlos- Section 5 presents the major results and 
analysis. A summary and conclusions are provided in Sec- 
tion 6. 



2 SAMPLE AND OBSERVATIONS 

Our s ample of stars was drawn from the IChristlieb et al] 
(2004) catalogue of FHB stars selected from the Ham- 
burg/ESO survey. As part of a large programme of follow-up 
observations of the stellar content of the HES and HK sur- 
veys, over 1100 FHB stars were observed at the UK Schmidt 
Telescope (UKST) 1 , using the six degree field (6dF) fibre- 
fed spectrograph. The 6dF instrument allows one to observe 
up to 150 objects simultaneously across the 6° field of the 
Schmidt telescope, obtaining spectra at 2 A resolution. More 
than forty nights of observations took place between April 
2001 and April 2003. Integration times were typically 4-5 
times 2700 s per field, depending on observing conditions. 

The targets were, in most cases, selected by HES field 
number, with roughly equal numbers of fibres used for metal- 
poor candidates and FHB stars. The data were reduced us- 
ing the Figaro and 6dfdr software packages, the latter being 
written specifically to handle the fibre data from the 6dF 
instrument. This process yielded FITS files containing the 
reduced spectra of all 150 fibres. These files were separated, 
with one spectrum per file and assigned a unique name based 
on the star name, field number, and date of observation. A 
hand-screening process was then used to exclude obvious 
problems and artefacts. 

Catalogues of FHB stars may contain a significant frac- 
tion of contaminants - higher surface gravity A-type stars 
which fall within the colour range of the survey (in the case 
of the HES, -0.2 < B—V < 0.3). Since the HES FHB candi- 
dates have been shown to have a much lower level of contam- 
inatio n from higher-gravity stars (< 16 %: Iciiristjiebetal, 
2004), as compared to the HK survey (~ 50%: lBeers et all 
1996), we considered only the HES stars in our observations. 
We made no attempt to separate the potential contaminants 
from the true FHB stars, judging that a ~ 16 % level is ac- 

1 The UKST is operated and supported by the Anglo- Australian 
Observatory, on behalf of the astronomical communities of Aus- 
tralia and the UK 
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Figure 3. Sky distribution of all 530 stars in 6 northern fields. 



Figure 4. Histogram of heliocentric distances for the final sample 
of 530 stars. 



ceptable for our present purpose. While it is, in principle, 
possible to further r efine the classif i cation , using the spec- 
troscopic method of IClewlev et all (|2002l ) would result in 
a sample clean to about the ~ 12 % level and hence noth- 
ing would be gained. Since systematic, accurate photometry 
is not available, methods based on this are not applicable. 
Conversely, we have no way of obtaining RVs from the ob- 
jective prism of the HES, but can measure velocities in the 
6df spectra, since the S/N requirements for RVs are less 
stringent than those for classification. 

2.1 Radial Velocities and Selection 

Sersic profile s were fit to the H7 and EM Balmer lines, as 
described bv IClewlev et af] J2002). Due to changes in the 
spectrograph settings over the long observing period, these 
were the only two prominent lines common to all spectra. 
We rejected any star for which we could not obtain a radial 
velocity measurement in both lines. We also noticed large 
systematic discrepancies (> 150 km s -1 ) between the mea- 
surements of RV from the H7 and H<5 lines in some stars. 
These stars were noted to be concentrated in one field (HES 
field 573), hence all stars in this field were rejected. One 
star was also rejected due to difficulties with sky subtrac- 
tion. This left a total of 827 stars. 

We sought to combine the velocity measurements from 
the two lines into an average value, but we first used the 
difference in the measurements of the two lines to provide 
an internal consistency check on the radial velocities. Fig.Q 
shows a plot of the difference between the H<5 and H7 line, 
versus signal-to-noise ratio. We applied cuts, retaining all 
stars with measurement difference less than 100 km s -1 and 
S/N greater than 15, as shown in Fig An average RV 
was then computed from all available measurements. He- 
liocentric corrections were calculated using the IRAF task 
rvcorrect. This selection yielded 636 RV measurements. 

Twenty-one stars in our sample were observed at more 
than a single epoch and thus enable a further check on the 
measurements. The correlation between the two average RVs 
for these stars is shown in Fig. [2] The difference between 
these two measurements for a given star - the range of mea- 
surements - provides an estimate of the standard deviation 
of the underlying error distribution <|Pearsonlll926l : iTippettl 



1925). This distribution will be normally distributed with 
\i — RV and a = error in the measure ment. Applying a cor- 
rection factor of 1.12838 (see Table X o flTippettll925i) to the 
absolute value of the range, we obtain a mean statistical ve- 
locity error of SOkms^ 1 . This agrees well with the assump- 
tion of Gaussian errors, under which we obtain 27kms _1 , 
dividing the standard deviation of the velocity differences 
by v2- We therefore adopt 30kms _1 as the formal error 
on our radial velocity measurements. To definitively charac- 
terise the error on a single RV measurement would require 
many repeated velocity measurements of the same star. For 
the analysis, the two RV measurements were combined into a 
single entry by averaging the two independent RVs, weighted 
by the S/N in the respective spectra, yielding RVs for 615 
unique stars. 



3 KINEMATIC MODELS 

The density distribution of the stellar halo follows a power- 
law, R~ a , with Galactocentric radius, R. The value of a 
has been the subject of many studies, using RR Lyrae stars, 
halo giants, and Blu e Horizontal Branch (BHB) stars (e.g. 
[ Preston et alJ Il99ll: iKinman et alJ Il994t ISluis fc Arnold! 
1998: iMorrison et al.l l200fj) . Mos t of these converge on a 
value of a « 3.0-3.5. In contrast. iRobin et alJ (|200Ch find a 
lower value consistent with the range a — 2.0-2.75, whilst 
noting significant degeneracy between a and the flatten- 
ing parameter c/a. This contrasts with the density distri- 
bution of the Galactic dark halo, which follows R~ 2 (i.e., 
the isothermal case expected from a flat rotation curve), at 
least in the outer parts where it dominates over the luminous 
matter. 

Stars in the halo of the Milky Way are drawn from 
a population with a total dispersion on the order of 
200 km s~ 1 , i.e. they have velocities needed to maintain equi- 
librium with the Galactic gravitational potential. In the stel- 
lar halo, this velocity is directed roughly equally into three 
components, resulting in a system with little bulk rotation 
and thus mainly "pressure support" against gravitational 
collapse. 

The three components of the stellar halo's velocity el- 
lipsoid have been measured accurately in the solar neigh- 
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Figure 5. Comparison of measured oxos *° * ne * wo models. The horizontal dashed line shows an isothermal model with the value 
measured by Sirko et al.. The solid line is the line-of-sight velocity dispersion predicted by the SLFC model for the average line-of-sight 
of the field. The error bars show the lcr errors. 



bourhood, for stars within roughly 1 kpc of the Sun, to be 
a = (<T R l<Te o 6 ) (140, 100, 100) kms -1 JChiba fc Beersl 
l2000l : iGouldl 120031 . where we have adopted a spherical co- 
ordinate system (R, 9, <f>) around the Galactic centre. This 
determination is based on stars for which all three compo- 
nents of the space velocity can be measured via both radial 
line-of-sight velocities and proper motions. Beyond a few 
kpc from the Sun, usually only the radial velocity can be 
measured (a situation which will chan ge dramatically when 
the ESA GAIA mission is completed: lPerrvmar] |2002) . and 
special techniques are required to reconstruct the three com- 
ponents of the velocity ellipsoid. One can either propose and 
fit data to models of the velocity ellips oid as a function of po- 
sition in the Galaxy - as performed bv lSommer-Larsen et alJ 
lll997l) for some ~ 70 stars - or attempt to recover the three 
components directly iSirko et al.ll2004a) . 

The velocity ellipsoid of the stellar halo is a consequence 
of its density profile and the total matter distribution in the 
Galaxy. Let us assume that the 3-D dispersion of the veloc- 
ity distribution is isothermal (i.e., it has the same magnitude 
in all three components), and that the density falls off like 
R~ a . It then follows, for a Galactic potential dominated by 
a dark halo with a flat rotation curve (characterised by con- 
stant circular velocity Vb) that a = Vo/^/a. For the Milky 
Way, Vb ~ 220 kms -1 , and thus a « 120 km s" 1 for each 



component teinnev fc TremainellT987l) . The total velocity 
dispersion is then of order ^f%a ~ 200 kms -1 . 

The stellar halo is clearly not isothermal in the vicin- 
ity of the Sun (the component of the velocity dispersion 
in the radial direction is « 40% larger than the two other 
components of the velocity ellipsoid); nevertheless, the total 
velocity dispersion is of order 200 kms -1 . 

What is t he velocity dispersi on like elsewhere in the 
halo? SLFC and Sir ko et ail (|2004a) have attempted a recon- 
struction from radial velocities of distant halo stars. SLFC 
proposed a model for the velocity dispersion cr(R, 9, <f)) as a 
function of Galactocentric radius R alone; the model per- 
mitted sharp deviations from isothermal behaviour (such as 
is seen in the solar neighbourhood) and was fit to both 
the local and distant halo stars. This model was subse- 
quently refined using a much larger sample o f abou t ~ 700 
distant halo stars by ISommer-Larsen et all lll997l) . It was 
found that the radial anisotropy (<tr > £>>,e) seen near 
the Sun persists into the inner halo (R < 8 kpc)), chang- 
ing to tangential anisotropy in the o uter halo (roughly 
R > 15 - 20 kpc). ISirko et al] J2004at) present an osten- 
sibly opposing view; from radial velocities of 1170 halo 
stars found in the SDSS, they found that the velocities 
of (distant) halo stars are very close indeed to isothermal 
(era tTfl.g-J = (1 01.4 ± 2.8, 97.7 ± 16.4, 107.4 ± 16.6) kms" 1 . 
The lSirko "etaU study is essentially of the outer halo, where 
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the SLFC model is not yet well constrained, hence the osten- 
sible contradiction with the non-isothermal model of SLFC 
is, as those authors argued, probably not significant. 

We treat the halo's kinematics statistically, charac- 
terised by the velocity ellipsoid alone; i.e. we explicitly ignore 
substructure, such as dissolving satellites or other features 
in density and/or velocity space. From the point of view of 
halo substructure, the halo can be divided into inner and 
outer regions. The inner halo is that part which is within 
~ 15 kpc radius from the Galactic centre. In this region, the 
relevant dynamical time scales are expected to be consid- 
erably shorter than the age of the Galaxy. This expecta- 
tion is borne out by the lack of detected substructure (e.g. 
iBrown et aDl2004h . Thus, in the inner halo, the stellar pop- 
ulation appears likely to be fairly well mixed, with a smooth 
density distribution. There is additional evidence that this 
is the case in the sol ar neighbourhood, based on an analysis 
of 4588 sub-dwarfs iGouldll20oA 1200411 . In the outer halo, 
beyond some 20 kpc or so, dynamical time scales approach 
the age of the Galaxy. In these regions it has long been sus- 
pected that the halo would contain the debris of dissolving 
satell i te galaxies; this is now well established (e.g. [Ma icwski 
Fool iThataet aUlTpl: iLvnden-Bell fc Lvnden-Belll ll <M 
iMorrison et alJ 2000). While the outer halo is undoubtedly 
quite lumpy, most of the stellar halo's mass, and most of 
that part probed via our data, is in the well mixed inner 
halo. Additional tests of this assertion will be the subject of 
a forthcoming study. 



4 FINAL SAMPLE AND VELOCITY 
DISPERSION 

We analysed the kinematics of our stars following the 
methodology of SLFC. This involved dividing the sample 
into regions on the sky in which there are sufficient stars to 
measure olos f° r t wo distance bins along the line of sight. 
The stars in each field were split evenly into the two dis- 
tance bins. The requirement we applied, that each bin should 
have at least twenty stars, forced us to exclude all stars with 
Galactic latitude b < 0° . This left us with a total of 530 stars 
with which to study the kinematics of the Galactic halo. The 
distribution on the sky of all 530 FHB stars in our final sam- 
ple is shown in Fig. [3] The average distance to the stars is 
6.7 kpc; the distance distribution is shown in Fig0] 

In calculating the velocity dispersion along the line-of- 
sight (olos) we excluded the upper and lower 5 % the data, 
in order to protect these numbers from the adverse effects 
of outliers. The result is adjusted by a factor of Q * 89 to 
recover the true standard deviation of the original d i stribu - 
tion, a process described in detail bv lMorrison et all |l990). 
The olos were corrected for a 30kms _1 velocity measure- 
ment error by subtraction in quadrature and then compared 
to the model predictions. We have consistently used veloci- 
ties in a heliocentric frame throughout the calculations, with 
the corresponding distance. Since the model predictions are 
made in a galactocentric frame, a projection factor was ap- 
plied to recover the appropriate line-of-sight value. 



Table 1. Definition of fields in Galactic co-ordinates, showing 
Bottom Left Corner (BLC), Top Right Corner (TRC), number of 
stars, distance and 0"los of stars in field. Errors are lo\ 



Field 


BLC 


TRC 


Num 


Dist 


CTLOS 




(l,b) 


(1,6) 




(kpc) 


km s — 1 


Field 1 


(220,45) 


(275,70) 


22 


6.3 ± 1.2 


108 ± 17 








22 


11.5 ±2.3 


110 ± 17 


Field 2 


(275,45) 


(320,70) 


78 


5.2 ± 1.0 


115 ± 10 








79 


9.0 ± 1.8 


125 ± 10 


Field 3 


(320,45) 


(360,70) 


31 


5.5 ±0.9 


92 ± 12 








32 


8.9 ± 1.7 


120 ± 15 


Field 4 


(220,20) 


(275,45) 


24 


5.0 ± 1.2 


104 ± 16 








24 


8.7 ± 1.0 


149 ± 22 


Field 5 


(275,20) 


(320,45) 


44 


3.9 ±0.8 


113 ± 13 








15 


7.5 ± 1.8 


120 ± 13 


Field 6 


(320,20) 


(360,45) 


64 


4.3 ±0.8 


84 ±8 








65 


8.2 ± 1.9 


98 ±9 



5 RESULTS AND ANALYSIS 

Table Q in addition to listing the fields, shows the observed 
velocity dispersions and the 1-a errors. These have been 
used to test both an isothermal, isotropic halo with veloc- 
ity dispersion of 1 01.6 km s -1 (i.e., the value observed by 
ISirko et al.l 12004^ ^1 and the predictions of the SLFC model. 
Fig shows this comparison for all six northern fields. The 
SLFC predictions are shown as the solid line, while the 
isothermal model is the horizontal dashed line. Visual in- 
spection leads us to conclude that both models fit the data 
equally well; a \ 2 analysis confirms this view. Both mod- 
els arise from predictions rather than fits to the data, and 
hence have no free parameters. For the twelve ctlos in the 
six fields we measured a reduced x 2 of 1.31 and 1.85 for the 
SLFC and isothermal models respectively, with 12 degrees 
of freedom each. Since this is not a statistically significant 
difference, both must be regarded as fitting the data equally 
well. 



5.1 Outer Halo Comparison 

Sirko etal ] J2004bl) have provided a similar sample of 1170 
FHB stars in the outer halo (distances up to ~ 100 kpc). 
From this sample they derive an estimate of the isothermal 
velocity dispersion of 101.6 kms -1 . We have analysed their 
data in the same manner as above (see their Table 3 for dis- 
tances and heliocentric red-shift). The results are shown in 
Fig ED At the distance of the Sirko data, both models pre- 
dict similar line-of-sight velocity dispersions. Therefore these 
data are not an effective test of the SLFC model. These data, 
divided into nine fields and three distance bins per field, fit 
both models equally well. The reduced \ 2 statistics for the 
two models are 2.37 and 2.19 for the SLFC and isothermal 
models respectively, for 35 degrees of freedom, excluding the 
obvious outlier in panel 3 (third from left, top row in FigEJ. 
This data point i s likely cause d by the Sag ittarius stream 
Jlbata et al.lll994h . as noted bv lSirko et all Again, the dif- 
ference in x 2 is n °t statistically significant. 
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Figure 6. Comparison of Sirko et al data oxos to the models fit by Sirko et al. (dashed lines) and SLFC (solid lines). Velocities 

are heliocentric and distances are in kpc. 



5.2 Discussion 

It is clear that neither this samp le of inner Galactic halo 
stars, nor the outer halo sample of lSirko et al] i200 4b1 . per- 
mit us to distinguish between the two competing models. It 
is worth noting h ere that the SLFC model is a physically 
realisable system (iFlvnn et al.lll996h . whereas the isother- 
mal model does not account for the locally observed halo 
velocity dispersion anisotropy. Where then should we look 
in order to best discriminate between these two models? 

We have used the model of SLFC to determine regions 
on the sky where the two models differ in their predictions 
of ctlos by more than 30kms _1 . This area is restricted to 
lines-of-sight toward the Galactic centre and at low Galactic 
latitude. Along these lines of sight one might expect a clear 
signature of the SLFC model, as demonstrated in Fig d A 
survey of 500 FHB stars at (I, b) = (340 to 20; -30 to + 30) 
with distance bins centred at 4, 8 and 13 kpc would provide 
a good discriminant between the models. This effect may 
also be seen in Field 6 of Fig |S] which is the closest the 
bulge. 



6 SUMMARY AND CONCLUSIONS 

Line of sight velocities have been measured for a sam- 
ple of 530 FHB stars in the inner Galactic halo, with 



an average heliocentric distance of 6.7 kpc. We have used 
the se velocities to test the veloc ity dispersions predicted by 
the lSommer-Larsen et ail lll994l) halo model, and compared 
with those predicted by an isothermal model. The former 
model allows for marked anisotropies to account for the well- 
constrained parameters of the local Galactic halo. These new 
data are equally well described by both models. We also 
compared the model predictions to the sample o f 1170 FHB 
stars in the outer Galactic halo, as published bv lSirko et al] 
(2004b). We confirm their findings that both models fit these 
data equally well. 

We suggest that the best place on the sky to differen- 
tiate these two models is towards the Galactic bulge and at 
low galactic latitudes; we propose that future efforts con- 
centrate on the more challenging task of recovering clean 
samples of halo FHB stars at low galactic latitudes. 
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